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For the zero temperature Glauber dynamics of the g-state Potts model, we calculate the exact 
distribution of domain sizes by mapping the problem on an exactly soluble one-species coagulation 



exact expression, we extract its behavior in various regimes. Our results are tested in a simulation 



model (^4+^4 — * A). In the long time limit, this distribution is universal and from its (complicated) 

o 

and compared to the predictions of a simple approximation proposed recently. Considering the 

: 

■ dynamics of domain walls as a reaction diffusion model A + A — > A with probability (q — 2)/(q— 1) 



and A + A — > with probability !/(<? — 1), we calculate the pair correlation function in the long 



■ time regime. 

o 
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1 Introduction 

When a ferromagnetic system is quenched from the high temperature phase to a temperature below 
its Curie temperature, one observes a pattern of growing domains. In the long time limit, when the 
typical size of domains becomes much larger than the lattice spacing (or the correlation length) but 
is still small compared to the system size, the domains form a (statistically) self similar structure 
(see [|| and references therein). For a non-conserved order parameter, it is well established that 
the size of the domains grows with time t like i 1//2 . Much less is known on the distribution of 
domain sizes. 

The purpose of the present paper is to give the exact distribution of domain sizes in the case of 
the one dimensional g-state Potts model evolving according to zero temperature Glauber dynamics 
(in Id, the low temperature phase reduces only to zero temperature). As [Q, |J] the average 
domain size grows with time like t 1//2 , the distribution of domain sizes can only be determined up 
to a change of scale and we will rescale the length x of the domains so that for the distribution 
g(x;q), the average domain size is unity (/ x g(x\q) dx = 1). So far, this distribution has been 
calculated only for q = oo [|], |(§ because it is related to the probability that two walkers do not 
meet up to time t. Our goal, here, is to extend this result to arbitrary q > 1. 

Our approach is an generalization of a calculation done recently |7], ||] to obtain the fraction of 
persistent spins (i.e. spins which never flip up to time t). We calculate the probability Q(N) that 
N given consecutive spins are parallel at time t and explain how g(x;q) can be extracted from 
Q(N). As for the number of persistent spins, the full expression of Q(N) is rather complicated 
but it can be used to write explicit formulae in various limits. 

It has been noticed for a long time that the zero temperature Glauber dynamics of the Id 
Potts model is fully equivalent to a single species reaction-diffusion model ||, ||, |(| ^, [K|. If one 
represents each domain wall by a particle A and if the initial spin configuration is random with 
no correlation, it is easy to show that the particles A diffuse along the line and that whenever 2 
particles sit on the same bond, they instantaneously react according to 



We will see that several of our results can be reinterpreted as properties of this Id reaction diffusion 




A with probability (q — 2)/{q — 1) 
with probability l/(q — 1) 



(1) 



problem. 



2 



An easy way to implement the zero temperature dynamics of the Id g-state Potts model is to 

say that during every infinitesimal time interval dt, each spin Si(t) is updated according to 

{Si(t) with probability 1 — 2dt 
S{-i(t) with probability dt 
Si+i (t) with probability dt 

This shows the close analogy with random walk problems [O]. As in M, [8|, |f2|, this analogy will 
be the basis of our calculation. 

The paper is organized as follows. In section 2, we recall several properties of random walks 
which are used in the following sections and the relation between random walks and the zero 
temperature Glauber dynamics of a spin chain. In section 3, we obtain, using several relations 
derived in ||, the probability Q(N) that, at time t, N given consecutive spins are parallel and we 
show how the distribution of domain sizes follows from the knowledge of the Q(N). In section 4, 
we give more explicit expressions in various limits and in section 5 we compare our predictions with 
the results of simulations. In section 6, we calculate the pair correlation for the reaction diffusion 
model (H|). 

2 Properties of random walks in one dimension and the q = oo 

case 

Let us first recall some properties of random walks in one dimension. Several of these properties 



are well known (in particular they can be obtained by representing the walkers as free fermions [13 



14j , |15| ) but we spend some time discussing them because they are essential to the understanding of 



the following sections. These properties will include the probability Cjj that two walkers starting 

(n) 

at positions Xi < xj do not meet up to time t, the probability c\ 2 2n that no pair meets up to 
time t between 2n walkers starting at positions x\ < 22 < ■■%2n- 

Consider a random walker on a Id lattice which hops, during each infinitesimal time interval 
to its right with probability dt, to its left with probability dt (and of course remains at the same 
position with probability 1 — 2dt). The probability pt(y,x) of finding the walker at position y at 
time t given that it was initially (at time 0) at position x evolves according to 

—pt(y,x) =Pt(y + l,x) +p t (y -l,x)- 2p t (y,x) =p t (y,x + 1) +p t (y,x- 1) - 2p t (y,x) (2) 
dt 



and the solution is 

m(ii. x) = 

2tt 



Pt(y,x) = I d9 cos(x - y)9 e-^ 1 —^)* . (3) 



Consider now 2 walkers on this Id lattice starting at positions Xj and Xj with x^ < Xj. The 
probability Ct(xi,Xj) = Cjj that these two walkers never meet up to time t evolves according to 

j t Ct{x, x) = C t {x + 1, x') + C t (x - 1, x') + C t (x, x' + l) + C t (x, x' - 1) - 4 C t (x, x') (4) 

with the boundary conditions that Ct(x,x) = at any time t and that Cq(x,x') = 1 if x < x'. It 
is easy to check that the following expression 

Ci j = Ct{Xi,Xj) = — / d6 J n — — e 4U COSO)t (5) 

J V ^ 27T./_ 7r 1-COS0 W 

satisfies (Q). It is remarkable that the probabilities of all the meeting events of N coalescing 
random walkers starting at positions x\ < X2 < •• < xn can be expressed in terms of the matrix 
Cij. For example, for 3 walkers starting at positions x\ < X2 < X3, the probability that none of 
them meets any of the other two up to time t is given by 

Cl,2 + c 2 ,3 - ci, 3 . 

For 4 walkers starting at x\ < X2 < X3 < X4, the probability that no pair meets up to time t is 
given by 

C l 2 2,3,4 = c l,2 C3,4 + Ci, 4 C 2 ,3 ~ Ci )3 C 2 ,4 • (6) 

Below, we will use the fact [|| that the probability that, up to time t, no pair meets among 2n 
walkers starting at positions x\ < x 2 < •••x 2n , is given by the Pfaffian of the 2n x 2n matrix Cjj 

C 1^2,...,2n-l,2n = ^"17 £ ( CJ ) C ^(l)^(2) " " " c a(2n-l),a(2n) ( 7 ) 

cr 

where the sum runs over all the permutations a of the indices {1, 2, • • • , 2n}, e(<r) is the signature 
of the permutation a and the matrix Cj.j is antisymmetrized 

Cjj = — c^j < when i > j . 

There are several ways of deriving (MR) , for example by considering the walkers as free fermions 



13, 14, 15, or by using the method of images. One can also simply write the equations similar to 
(Q) which govern the evolution of these probabilities and check that (||,|7]) do satisfy these equations 
with the right boundary conditions when Cij is solution of (0). (note that the above relations would 
not be valid in dimension higher than 1: it is only in Id that whenever 1 and 2 do not meet, it 
implies that 1 and 3 do not meet because of the order of the walkers along the line). 



Apart the probabilities that the walkers never meet, the knowledge of the matrix Cj allows 
one to calculate all the meeting probabilities between coalescing random walkers (8|. For example, 
when four walkers start at positions x\ < X2 < X3 < X4, one finds for the probability "012,34 that 1 
and 2 coalesce, 3 and 4 coalesce but 2 and 3 do not coalesce before time t 

(2) 

■012,34 = Cl, 3 + C 2 ,4 + 4,2,3,4 ~ c l,2 ~ C2,3 ~ C 3 ,4 (8) 

and for the probability "01,23,4 + ^1,2,3,4 that 1, 2 do not meet and 3, 4 do not meet (without saying 
whether 2 and 3 meet or not) 

(2) 

"01,23,4 + ^1,2,3,4 = Cl j3 + C 2 ,4 + 4,2,3,4 ~ c l,4 ~ C2,3 • (9) 

In the long time limit, all the above expressions (f],|7],||,||) remain valid with aj replaced by the 
asymptotic expression of 

c i j = Ct(x i ,x j ) = f(%^) (10) 

f{z)=^= r e -« 2 du. (ii) 

V 71 " •'0 

2.1 Relation between random walks and the dynamics of the spin chain: 

As in [l^, [7|, g|, we are going to use the close analogy between the problem of coalescing random 
walkers and the properties of a spin chain evolving according to zero temperature dynamics: 
assume that we want to calculate, for a Potts chain with no correlation in the initial condition, 
the probability that spins at position x% < X2 < ■■■ < xn are equal at time t. One can consider 
coalescing random walkers starting at positions x% < X2 < ■■■ < xn and if P(m,t) denotes the 
probability that at time t, there are m walkers left in the system, one has [[?], ^] 

N 1 

Probability^ (t) = S X2 (t) = ...S XN (t)} = £ Pfa,*)-^ . (12) 

m=l ^ 

This expression can be understood by noticing that initially the spins are random and that the 
probability that m spins have the same color in the initial condition is just l/q m ~ l . 

2.2 The case q = oo: 

In the limit q — > oo only the first term of the sum contributes and one finds that 

Probability{S xl (t) = S X2 (t) = ...S XN (t)} = 1 - C t (x u x N ) . (13) 
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This expression can easily be understood: for q = oo, all the spins in the initial configuration are 
different; if S xi (t) = S XN (t), this means that the two walkers starting at x\ and x^ have met 
before time t and therefore, at time t, all the spins located between x\ and xn are equal. It can 
be used to calculate the distribution of domain sizes in the q — > oo limit: if pi (Z) is the density of 
domains of length I, one has 

Probability^, (t) = S X2 (t) = ...S XN (t)} = £ (l-x N + Xl ) Pl (l) . (14) 

Z>X]\r— x\ 

This leads to 

pi(i) = 2C t (0, 1) - C t {0, l + l)- C t (0, 1 - 1) (15) 

and in the long time limit 



8r \y/8t 

For q = oo, this allows one to calculate the average domain size at time t 

and the distribution (/(a;; oo) of domain sizes || (normalized such that the average size is 1) 

For similar reasons, it is easy to see that for q = oo, the probability that S Xl (t) = S X2 (t) ^ S X3 (t) = 
S X4 (t) is given by 

Probability^ (*) = S X2 (t) ± S X3 (t) = S X4 (t)} = V>i 2 , 3 4 (19) 

This can be used to obtain ^2(^1^2)) the density of domains of length l\ followed by a domain of 
length I2 (because for q = 00, whenever two spins have a certain color, all the spins between them 
have the same color) 

Probability^ (t) = S X2 (t) ± S X2+1 (t) = S X4 (t)} = £ £ V2(hM) ■ (20) 



and consequently (|-11), the (normalized) distribution g2(x,y;oo) of neighboring domain sizes is 
given by 

92(x, y; 00) = (if = l( x + y ) [ e -(^W* _ e -(^/4l (21) 
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This means that the lengths of consecutive domains in the limit q —* oo are correlated. In partic- 
ular, if l\ and I2 are the lengths of two consecutive domains, one has 

(ilia) f°° , f°° , , v 3 



(0 s 



/ dx I dy g 2 (x, y; 00) x y = - ^ 1 (22) 



Remark 1: The knowledge of the distribution of domain sizes and of all the correlations between 
consecutive domain sizes for q = 00 gives the distribution of domain sizes for arbitrary q. If 
h,l2, ■■■■l n , ■■■ are the lengths of consecutive domains when q = 00, one can obtain the length I of 
a typical domain for finite q by 

1% with probability 

l\ + I2 with probability 

9 ( 23 ) 

h + h + ■■■In with probability 

This is because one way of performing the dynamics at finite q is to first make the system evolve 
with all the spins different in the initial condition (as in the case q = 00) and then suppress the 
domain walls with probability 1/q. This is due to the fact that for finite q, two initial values are 
identical with probability 1/q. 

A consequence of (^) is that, in the long time limit, the average length for finite q is given by 
(0) 

(l)g = -^(Ooc - -^—V2ri (24) 



It is clear also from (23) that the calculation of any other moment (P) of the length would require 
the knowledge of all the correlations between the l{. 

Remark 2: If these correlations between successive domain sizes for q = 00 were absent, the cal- 
culation of the normalized distribution g(x; q) for arbitrary q would be straightforward: neglecting 
the correlations in ( f23| ) would give for the generating function of g(x; q) 







dxg{x; q )e^l^dx = (Q-^Sfdxgjx^dx 

q — J dxg{x; oo)e ax dx 



(25) 



Then using the expression (18), one would find (by analyzing the limit a — > —00) 



/ \ 71 q TT 2 (3q-l)q 2 3 vr 3 (15q 2 - 6q + l)q 3 5 - 

q(x:q) = x ; r-5— x H ; r= x + U(x ) (26) 

yy 2q-l 24 (q-l) 3 960 (q - l) 5 V ' V ' 

(and by analyzing the pole in a) 

g(x; q) ~ exp [—A(q)x + B(q)] for large x (27) 



where A{q) is the root of 

irx 2 (q — l)x 



2 Jo 

and B(q) is given by 



/■oo 

/ dx xexp 



+ ^A(q) 



(28) 



q 2 irA(q) 
q7T + 2(q-l)A(qf 



m=**{ „,ir™ A ,~ • (29) 



For q = 2, this is precisely the prediction recently given in |16|, based on the assumption that the 
intervals l±,l2, • • are uncorrelated. We have already seen (21,22) that this assumption is not valid. 



Our goal in the following sections is to obtain the true distribution g{x; q) where these correlations 
have been taken into account. We will see that the exact expression of g(x; q) disagrees with 
(HJ2|J1|. 



3 Arbitrary correlation functions 

Our solution for the distribution of size of domains is based on our ability to write exact expressions 
for all correlation functions (valid at an arbitrary time) when the initial condition is random with 
no correlation (i.e. each spin takes initially one of the q possible colors with equal probability). 
These exact expressions involve the probabilities given by (||]6],@). One can show in particular 
that the probability that N spins located at positions x\ < X2 < ■■Xn are identical at time t is 
given by 

N-i 

Piob{S xl (t) = ...S XN {t)} = l-nYl ci,<+i + M 2 Z)^5+ijj+i - ^ J2 %i+i,j,j+i,k,k+i + ■■■ 

i=l i<j i<j<k 



"A I MC1.7V - V 2 C l},i+1,N + V 3 J2 C lli+l,j,j+l,N \ ( 30 ) 

[ i i<i 

where 



\ = q-l (31) 

M=^. (32) 
q 2 

(note that as long as iV is finite, ( |30| ) is the sum of a finite number of terms and is therefore a 
polynomial in the variable n). 

The proof of ( p0| ) is exactly the same as the one given in |7j : one shows that fl30| ) is equivalent 
to (|l^) for N coalescing random walkers starting at positions x\ < X2 < ■■■xj^ by considering all 
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the possible coalescing events between the N walkers. The weight a m in ( |30| ) of an event where at 
time t the N walkers have merged into m walkers is given by 



(m — 2)(m — 3) 9 (m — 3)(m — 4)(m — 5) ■> 
a m = 1 - (m - 1)// + ^ V - ol V 



2! 



3! 



. 9 (m — 4)(m — 5) q 
-A ( /i - (m - 3)^ 2 + i ^ -// 



(33) 



which is just a m = l/q m (see (plf ) and (p2[) ). This therefore proves that ( p0|) is equivalent to 



In ||, several ways of rewriting ( |30[ ) were given. In particular it was shown that ( |30| ) can be 
rewritten as the square root of a determinant 



Probst) = .. 



/ 1 \ 
10 
. 



V 



( 10 A \ 
-10 10 
-10 1 



1 / 





V -A 



1 

■1 / 



V 



By choosing N consecutive sites for x\ < X2 < ..xn, the expressions (30,|3_J) give the probability 
Q(N) that N consecutive spins are parallel and the density pi(l) of domains of length I follows 
from a relation similar to (|i~4| ) 

Q(N) = J2(l-N+l) Pl (l) . 

1>N 



The long time regime 

In the long time limit, Cj.j varies slowly with the positions x% and Xj of the sites, and becomes a 
continuous function c(xi,Xj) of these two positions. Moreover, for large N, all the sums in the 
expression ([}(]) become integrals. As in ||, one can show that in this continuum limit, (^) or 



(34) can be rewritten as 



Q(N) 



yJl-Hc(N,N)-\yJ-»c(N,N) 







exp 


-tr log M 

.2 



where this time c{x,y) = f{Vj=) as in ( ^0|Jll| ), the matrix M is defined by 



(35) 



M(x, y) = 6(x -y) + 2fi —c(x, y) 

ax 



and the quantities appearing in ( p5| ) are given by 



tr log M 



£ 

n=l 



(-2/u) 



n r N 



n 



dx 



l • • • 



N d d 
dx n —c(x 1 ,x 2 ) ■ ■ ■ —c(x n ,xi) 



(36) 



and 



c(N,N)= dyc(N,y)M- 1 (y,N) = y2(-2 f i) n dx x --- dx n c{N, Xl ) — c(x 1} x 2 ) ■ ■ ■ — c(x T , 
Jo ~ x Jo Jo ax dx 

(37) 

If one normalizes all the distances on the lattice to make the average domain size become 1, 
one defines x by 



N = V2iTt 



1 



and the distribution of domain sizes g(x;q) is then given by 

Q(N) = / dy(y-x)g(y;q) 



so that in the long time limit, one ends up with 



dy(y - x)g(y;q) 



1 — fiAi(x) — Ay— iaA\(x 



where 



and 



with 



and 



A i( z ) = ( _2 A i ) n / dxi--- dx n j(z,xx) — j(x 1 ,x 2 )---— j(x n ,z) 
^ Jo Jo dx 



71=1 



dx 



. . , l ^ {-2fi) n r z , r z , d . . d 

A 2{z) = --y / dxi--- dx n — 7(xi,x 2 ) • • • — 7(x n ,xi) 

9 ^— ' n Jo Jo dx 



71=1 



dx 



7(x, y) = F(y - x) = —= I e u du 

V 7 *" Jo 



A 7 (x,y) = -^ I exp 



-7r(y - x) V 



4(g-l) 2 

Differentiating (|39|) twice with respect to a; leads to the exact expression of g(x; q). 



(38) 



(39) 



(40) 



(41) 



(42) 



4 Expansions for the distribution of domain sizes 



The exact expression (^) for g(x; q) requires the calculation of A\ and A 2 which we did not 
succeed to do explicitly for general q and x. However, it can be used to expand to arbitrary orders 
in various limits: small x, large x, large q. 
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4.1 small x 

The small x expansion is straightforward. One just need to expand j(x, y) in (fi"^) in powers of 
y-x 

7 (x,y) = F(y - x) = ^(y -x) - £ [-^\y - xf + ^ (^) 5 (y - *)*... 



and to replace 7 by its expansion in ( j39| , [4C| ,41). Mathematica performs this to high order, and we 
found for g(x; q) up to order x 13 

_ it q 7T 2 q 3 3 tt 2 q 3 4 tt 3 g 5 5 

" 2g-l X 8 (g-l) 3 ^ + 24( 9 -l)4 X + 64(g-l)5 X 
^ g 5 6 tt 4 g 7 7 tt 4 g 6 (3 + 23g) 8 
120(g-l) 6X 768 (q-l) 7X 26880 (g - l) 8 X 

tt 5 g 9 x9 tt 5 g 8 (ll + 29g) 1Q tt 6 g 11 n 
12288 (g-l) 9X 483840 (g - l) 10 X 245760 (g - 1) 11 °° 
tt 6 g 10 (429 + 547g) 12 
+ 170311680 (g - l) 12 X 

tt 6 g 10 (64- 64g + 3157rg 3 ) 13 
+ 1857945600 (g - l) 13 °° 

+ 0{x lA ) . 

This expression is clearly different from the one obtained (26) in the uncorrelated approximation 



16]. Note that only the first term in both expansions is identical. 



4.2 Large x 

For domains much larger than the average size, we expect that, for finite g, g(x; g) decays expo- 
nentially 

g(x; g) ~ exp[-A(q)x + B(q)\ . (43) 



This can be understood by recalling that a large domain of size L at finite q is created (23) by 
combining n domains at q = 00 with n typically L/(l) q=00 . The probability of this is exponentially 
small in n. If the domains at q = 00 were not correlated, the constants A{q) and B(q) would 
be given by ( |28]p9| ). The existence of correlations between the lengths of domains ( p2[ ) alters the 
coefficients of the decay, but not the exponential decay itself. 



The large x behavior of g(x; q) is directly connected (39) to the large N behavior of Q(N) 



POO 

Q{N) = dy(y- x)g{y; q) ~ exp[-^(g)x + B(q) - 2 log A{q)) . (44) 

Jx 
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Therefore, we need the large x behavior of A\(x) and ^(a;) given by (40,41). It turns out that 
(as in ||) the prefactor >/l — fj,A\(x) — \y/—LiA\{x) vanishes as x — ► oo when q > 2 but has a non 
zero finite limit when q < 2. One has therefore to examine these two cases separately. 

• For q < 2, the prefactor (which has been calculated in equation (35) of ||) has the following 
limit 

yjl- nAi{x) - X^J-fiA 1 (x) -> \Jq(2 - q) for x -> oo . (45) 
To evaluate the large x behavior of A2(x), one can use well known results on Toeplitz deter- 



minants [17, [HI, namely that for an even function a(u) which decays as u —* =too 



7^ / dtti • • • / (itin a(u\ — U2) ■ ■ ■ a(u n — u\) c± 

„=i n ^ 



— / log[l + 2fia(k)]dk + / u du—^ 
2vr j_oo Jo 47r z 



e~ iku log[l + 2fid(k)]dk 



(46) 



where 



e iku a{u)du 



When this is used in the case of A2{x) given by ( pi]) and ([42]), one finds that 



a(u) 



g-1 6XP [4( g -l) 2 
( -k 2 (q- 1) 2 1 



-7ru 2 q 2 



and 



a(fe) = — 2exp 



7rg z 



(47) 



This leads to the following expressions in terms of fx — — - 

1 q 



A(q) 



r 

(4/i)« 



4 q 



1 ^ jr ! 2 



n=l 



B(g) - 2 log A(q) = - log g(2 - q) + — £ i^- £ 
2 47r ^— * n 



„_ 2 p=1 ^p(n-p) 



(48) 



(49) 



which can be rewritten in order to make each term have a finite limit in the previous expres- 
sion when q — » 2 



B(g) -21ogA(g) = logg 



g - 1 1 ~ (4/i)" 



n-l 



1 



g 2 47r 



n=2 



71 



-7T + - 

p= i VP(n - p) 



(50) 



In ( f48|) and (|5^), the series are convergent for all values of q (since (q — l)/q 2 = fi < 1/4) 
and so these expressions can be used to calculate numerically A(q) and B(q). In section 5, 
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we will compare the values predicted by (|4q) and (|5(]) to those ( ^q , 2£ ) of the approximation 
where the domains are uncorrelated. 



Remark 1: the radius of convergence of the sums which appear in (48,^) is /U = 1/4. 
So as /i — > 1/4, that is g — > 2, these expressions could become singular. As [i < 1/4 when 
q > 2, ( f48[) and ( |50|) can be computed for both for q < 2 and q > 2. However, it is easy 
to see that these expressions as a function of q are not analytic at q = 2 (for example it is 
possible to show that -^^-A(q) —* y/n/i as (? — ► 2 so that ^4(g) given by fl4£j) has a cusp at 
(7 = 2). We will argue below, that for q > 2, expressions ( |4S| ) and (|5^) are no longer valid as 
they are but should be replaced by their analytic continuation from the range q < 2. 

Remark 2: a priori, in the Potts model, q is an integer and so the case q < 2 is of lit- 
tle interest. However, Fortuin and Kasteleyn [20] have shown a long time ago that non 



integer values of q have physical realizations as cluster models. Here, it is very easy to check, 
that if one considers the zero temperature dynamics of an Ising chain where the spins are 
initially uncorrelated but a non zero magnetization m, the distribution of the sizes of domains 
of + spins is exactly the same as for the q state Potts model when 

2 



1 + m 

In particular, q — ► oo corresponds to an initial condition where the + spins are very rare 
whereas q — ► 1 corresponds to very few — spins. So m can take any value between —1 and 1 
and q can vary continuously between 1 and oo. 

Also, for q > 2, the reaction diffusion model (||) makes sense for a non-integer q. 

For q > 2, the large x behavior of A<i{x) is exactly the same as for q < 2. However the 
prefactor in ( |39| ) vanishes as x — > oo. The evaluation of this prefactor for large x is not 
simple. A similar calculation was done in || and the result for q > 2 turned out to be simply 
the analytic continuation of the result obtained for q < 2. Here we assume that this property 
remains true for g(x; q) (in the small x expansion, at least, one can easily check that all the 
coefficients can be analytically continued at q = 2) and that A(q) and B(q) in the range 



q > 2 are given by the analytic continuation A(q), B(q) of their expressions (|48| , |50D for q < 2. 
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In the appendix, we give a way of deriving the following expressions of A(q) and B(q) 



o — l 4 a — 1 i 



4 o - 1 ^ n 3 / 2 

3 n=l 



(51) 



1 1 00 (4 ) n ri— ^ i 
Sfo) - 2 log ]%) = - log q(q - 2) - log 4 - log(- log 4 M ) + — ]T i-^- ]T - 

2 4vr ^ n ^ ^/p(n - p) 



1 



n=1 Tly 7T J^f-nlogin 



dve~ v .(52) 



4.3 The large q expansion 



In the large q limit, fi, A and 7(2;, y) can be expanded in powers of 1/g (see (31,32|f42|)), This leads 
to 

A 1 (z) = -^7 2 (0,z)+4 M 2 7 (0,z) £dx<y(0,x) + O ( 4 

A 2 (z) = -- + ofl 



and to 



Q(JV) = 1 - / due-™ 2/4 



+ T due-™ 2 / 4 ) + 2 f due"™ 2 / 4 - xe~™ 2 / 4 - 2 f dy [ V dze-^ z2+ ^ x ~ y ^} 

q I Jo Jo Jo Jo J 



+0 I T2 I • 



(53) 



which becomes using (3S) 



2 2 2 



£dze-X*-* ) 2 +* 2 W 4 j+O (^L 



(54) 



This gives for the moments of the distribution of the distribution g(x;q) 



{x n } 



2_\ n 

7T 



, n 



^ - ( " +i)r e +i )-^(T-)^5(;) r (^) r m 



1 



(55) 



and in particular, 
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Remark: because of (23) all the results to order 1/q can be recovered from the knowledge (18,21) 



of g(x; oo) and 52 (%, x'; 00) and it is easy to check that 

1\ /. 1 \ 1 f( 1 +^> , / 1. 

1 I 9 (1 + ~) x 'i 00 + - / ds 52 (1 + — )ac - z, z; 00 

qj \ q J q Jo \ q 

is equivalent to (|54|). It is straightforward in principle to generate higher order terms in ~ but the 



g{x;q) = ( 1 + - 



expressions become quite complicated. 

5 Simulations 

In order to check the validity of our results against the predictions of the independent interval 



approximation 16], we made a simulation. Because of (j23), the distribution of domain sizes for all 
values of q can be extracted from a single simulation done at q = 00. One starts with a random 
initial condition at q = 00 (this is done by choosing initially each spin with a different value: for 
example spin i is given color i) and let the system coarsen according to zero temperature Glauber 
dynamics. Once the system at q = 00 has evolved for a certain time, one obtains a system for an 
arbitrary value (integer or non-integer) of q by removing each domain wall with probability 1/q. 

The easiest properties one can measure are the moments of the distribution of domain sizes. 
Our data (for a system of 10 6 spins) for the second moment (x 2 ) of the distribution g(x;q) of 
domain sizes for q = 2, 3 and 5 are shown in figure 1 as a function of the average size (I) 00 of the 
domains for q = 00. For Glauber dynamics, this average (l)oo increases with time like t 1 / 2 . The 
larger {l)oo is, the closer we are from the asymptotic regime. However, as the simulation is done 
for a finite system, the results become noisy if one waits too long simply because the number of 
domains is too small to allow good statistics. 

These numerical results can be compared to the approximation where the correlations are 
neglected ( |25| ) 

( x ) independent = ~ ~ + ~ ■ (^) 
We see in figure 1 that our numerical data seem stable and accurate enough to show a clear 
discrepancy with the approximation. Unfortunately, our exact expression (|39| ) is complicated and 
we were unable to obtain closed expressions of the moments for the whole range of q. We could 
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test however our conviction that correlations between the domains at q = oo have an effect. If in 
the calculation of (x 2 ), we neglect the correlations between all pairs of domains except for nearest 
neighbor domains (|22|), we find an improved estimate 

/ 2\ _ g-14 2 q-1 /3 \ 

(x )i mprovec i - + ~ + 2— ^2~~ ^~ ~ 1 J ( 57 ) 

which, although not exact, seems to be in much better agreement with the results of the simulations. 

As we could not calculate the moments from (|39|) , we tried to test our results by measuring 
the exponential tail of the distribution g(x;q). It is always difficult to measure numerically a 
probability distribution because one has to use bins: if they are too narrow, the data are noisy 
and if they are too broad, all the details of the distribution are smoothened out. To avoid this 
difficulty, we measured the integrated distribution h(x; q) 

roc 



and we compared it to our prediction (|43| ) for the tail (48j50|) for q < 2 and (51,52) for q > 2. 



We also compared the results of our simulations to the prediction of the independent domains 



approximation (28,^9). In figure 2, we plot the logarithm of the integrated distribution divided 



by its predicted form [A(q)x — B(q) + log A(q)] + log h(x;q) versus x for q = 5 for a system of 2 



millions spins, after 1000 updatings per spin. We see that the agreement is very good when (51,52) 



are used whereas for (g8L 29 ) , there is a small but visible discrepancy. Of course, for x too large, 
our data are too noisy and for x too small, the asymptotic form has no reason to be valid , so that 
the range where the agreement is the best is 1 < x < 3.5. 



In table 1, the exact values of A(q) and B(q) obtained from (|48|j50| ) for q < 2 and (|5l| , p2|) for 



q > 2 are compared to those coming from the independent domains approximation ( 28| and |29|) . 
We see that, although different, they are very close. 

6 The reaction diffusion model 

The relation between the spin dynamics and the reaction diffusion model (||) can be used to 
calculate all kinds of properties of the reaction diffusion model. For the zero temperature Glauber 
dynamics of the Potts chain, it results from the analogy to random walks (section 2) that 

Probability^^) + S X2 (t)} = (— ) d,2 
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and that @ for x\ < %2 < £3 < x 4 
Probability {S Xl (t) / S X2 (t) and S X3 (t) / 5^ 4 (t)} 



2 r 



Cl,3 + C2,4 - C 2 ,3 - Ci )4 + c\ 



(2) 



2,3,4 



In the long time regime, this becomes using (pX|) and ( |TTD for x\ = 0, 22 = 1, £3 = r, 24 = r + 1 



Probability{S (t) / Si(t)} ~ #'(0) 



Probability{5 (t) / Si(t) and 5 r (t) 7^ S r+ i(t)} ~ 



9-1 



-#"(r) + H"(r)H(r) - (H'(r)) 2 + (#'(0)) 2 



where 



2 r/v's* 2 

~ ' e~ u du 



2 r 

V7T JO 

If these expressions are interpreted in terms of domain walls (i.e. of particles in the reaction 
diffusion model), one obtains that the density p of particles and the probability ^(O, r) of finding 
a pair of particles at and at r are given by 



-H'(0) 



i2id 



p 2 (0,r) 



1 



-H"{r) + H"{r)H{r) - {H'{r)f + (H'(0)) 2 



By eliminating the time t between these two equations one finds that in the long time limit, 
the probability of finding a particle at and r is given by 



P2(0,r) =p 2 



-2z* , o„„~* 2 / e ~u* du 



l-e~ 2z +2ze 



(58) 



and 



q 

Z = r— — P T 



q-1 2 

We see in particular that the correlations decay like a Gaussian instead of the "usual" exponential. 



25] that several properties of reaction 



Remark: it has already been noticed [21, 22, 
diffusion models can be calculated exactly, by writing closed equations for the probability that an 
interval of length I is empty (or two intervals are empty...). This method (which is rather close in 
spirit to our approach as it essentially considers that walkers which do not meet in one dimension 
are free fermions) could also be used to recover 
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7 Conclusion 



The long time regime of domain growth phenomena is by many respects analogous to what happens 
near a critical point in a second order phase transition: this regime is characterized both by 
universal exponents and by universal scaling functions. In the present paper, we have determined 
the distribution of domain sizes in this regime, for the (/-state Potts model in one dimension. This 
distribution is universal, in the sense that, at least, it would remain the same for short ranged 
correlations in the initial condition. 

The exact expression we found is rather complicated (|39|^) but it can be simplified in several 
limits (section 4). Our exact results are different (although rather close to) the predictions of a 
simple approximation proposed recently [16]. However, with sufficient numerical effort, we think 
that our simulations of section 5 indicate the validity of our results against those (^£-29) of that 
approximation. 

The present work has a lot in common with ||] where the exponent characteristic of the per- 
sistent spins have been calculated exactly for the same system. Here we are concerned with space 
properties whereas ]|] was devoted to time properties. The existence of non trivial exponents for 



the number of persistent spins 27, 28, 29, 30] in d > 1 has certainly its counterpart in terms 
of distribution of domain sizes. However, in d > 1, the size is one among many ways of character- 
izing a domain (its volume, its perimeter, its shape |j31j| , its number of neighboring domains...). It 
would be interesting to know whether the probability distributions of all these characteristics of 
domains are universal in d > 1. Another open question for systems in d > 1 is how to define the 
domains in presence of thermal noise (i.e. at non-zero temperature) and how to distinguish them 
from thermal fluctuations. When the typical size of domains is much larger than the equilibrium 
correlation length, it seems intuitively easy to make this distinction, but to our knowledge, a clean 
way of measuring the size of domains (in presence of noise) has not yet been proposed. 



We thank Vincent Hakim and Vincent Pasquier for useful discussions. 
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8 Appendix : Analytic Continuations 



In this appendix we show how to analytically continue the expressions (48, [5C]) of A{q) and B(q) 
obtained in the large x expansion for q < 2 to the range q > 2. It is easier to start with the 
expressions coming from ([44|-f4"6|) , that is 



A{q) 



Air 



dk log[l + 2/j,a(k)} 



1 l r°° 

B(q)-2\ogA(q) = -\ogq(2-q) + — i y dy 

2 87T Z Jo 



dke iky log[l + 2/m(fc)] 



where d(k) is given by (|47|) 



a{k) 



-2 exp 



-k 2 {q-lf 
irq 2 



(Al) 
(A2) 

(A3) 



As q — > 2, that is [i — > 1/4, the integrands in ( |A1| ) and (A2) become singular at k = 0: two complex 
zeroes in the complex plane of the variable k approach k = as q — > 2 and the exchange of these 
two zeroes is at the origin of the analytic continuation. 



A generic example: to make the analytic continuation of ( |A1| ) and ( |A2| ) simpler to understand 
we treat first where the expression ( |A3| ) of d(k) would be replaced by 

-2 



d(k) 



(A4) 



l + Q(k) 

where Q(k) is a real polynomial of degree 2n with (Q(k) ~ k 2 as k — > and such that Q(k) + 1 
and Q(k) + 1 — 4^ have both only complex zeroes). One can therefore write 

4/x 



1 + 2/jd(k) = 1 



n (Zj .- fc )( z *_ fc ) 

l + Q(k) fj^yj-k^-k) 



(A5) 



where the yj are the n zeroes of 1 + Q{k) with positive imaginary parts and the Zj are the n zeroes 



of 1 — 4/i + Q(k) with positive imaginary parts. One can rewrite (A1,A2) as 

1 r°° „ 2u k a'(k) 1 Z 100 „ An k Q'(k) 



. 00 „ 2fi k a (k) 

A(q) = — / dk— Vtt - — 

w 4vri_ 00 l + 2^a(A: 4vr 



1 1 r°° dn 

B{q)-2\ogA{q) = -\ogq(2-q) + — 1 / J- 

2 8tt z Jo y 



dk 



[l + Q(fc)][l-4/z + Q(fc)] 



dk e iky 



AfiQ'(k) 



[l + Q(k)][l-4fi + Q(k)} 



(A6) 



(A7) 
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and find using the theorem of residues 



3=1 



and 



B(q) -2logA(q) = - log q(2 



j=ij'=i 



IZ; 



iz-, - iyj 
IZ*, - IZj ^ 



(A8) 



(A9) 



These expressions (A8) and (|A^) are valid for q < 2. To obtain their analytic continuation 



A(q),B(q) to the range q > 2, we have essentially to exchange the role of z\ and z*, the two zeroes 
of 1 — 4/i + Q(k) which go through k = when q = 2. So for q > 2, 



(A10) 



and 




(All) 

where we have added an infinitesimal imaginary part to q to define the analytic continuation of 
log(izi — izl) and log q{2 — q). 
Using the fact that 

n 

!]>! - yj )( Zl - y*) = 1 + Q( Zl ) = 4/i 

3=1 



and that 



Y[( Zl - z j )(z 1 - z*) 

3=2 



Q'(zi) 

z\ - z\ 



dH) b 



ecomes 



B(q) -2logA{q) = -logq{q-2) + log4/i - log |Q'(zi)| + log|*i - z\\ 



N 



■ £ log {(zt - yj ){ Zl - y*)} + J2 log {(zi - z*)(zt - Zj )} 

3=1 3=2 



1 



i=ij'=i 



ly 



3 1 



IZj 



iZj, - iyj 
iz*, - IZj _ 



This can be further transformed using the following identity 
dk 



2-Ki 



z\ 



1 \, ( (k - zj)(k - z*) 



(k- yj )(k-y*) 



log 



(A12) 



'{z\ - z*)(z* - ZjY 
(zi -y*){z{ 
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to become 



B(q) - 2 log A(q) = - log q(q - 2) + log4/i - log |Q'(^)| - log \ Zl - z\\ 



+ 2-ni Xoo ^ (& — Z\ k — z\ ' ' 



Lastly, using the identity 



one finds for q > 2 



4/i 



2/xa'(zi) 



1 + Q(fc) 

IvviogJf^-"^ 

W(^i) _ Q'(z x ) 



(A13) 



5(g) - 21ogA(g) = -logg(g - 2) - log 1 2^(zi)| - log|*i - sj| 



1 



1 



1 f 00 

+ / dk . 

2m J-oo \k — z\ k — z\ 



log (1 + 2fid{k)) 



+ ib[ ydy 



dke iky log[l + 2fm(k)] 



This expression together with ( |A10 ) 

i i r<x> 

A(q) = ~{z{ -zx)-— j ^ dk log[l + 2/xa(fc)] 



(A14) 



(A15) 



give for the generic example ( |A3| ) the analytic continuation A(q),B(q) of A{q) and -B(g) to the 
range q > 2. 

Although (|Al4l|Al5|) were derived for 5(fc) of the form ( |A3[ ), it should remain valid for a broader 
class of functions, probably as long as d{k) is such that only a pair of zeroes of l + 2/ia(fc) exchange 
as q crosses 2. 

We believe that flAliD and flAlgp remain valid for d{k) given by (4^) so that 



• 9 \/tt / — -, -. — 

z\ = i -V- log 4/i 

g- 1 

g — 1 4 g — 1 i n 6 l z 

* * n=l 



(A16) 



B{q)-2 log A(g) = - log g(g - 2) - log 4 - log( 



n=2 



«=i \W - p) 



1 



dv e~ v (A17) 
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Figure caption. 

Fig. 1. The second moment (x 2 ) versus (l)oc for q = 2, q = 3, q = 5 for a system of 10 6 spins. 
The dashed lines indicate the values predicted by the improved approximation ( |57|) and the 
plain lines those of the independent approximation (|56| ) . 

Fig. 2. The logarithm of the (measured) integrated distribution h{x;q) divided by its predicted 
asymptotic form A(q)x — B(q) + log A(q) + log h(x; q) for a system of 2 x 10 6 spins after 1000 
updatings per spin. When the exact expressions of A(q) and B(q) are used, the agreement 
is very satisfactory, whereas when on takes the values predicted by the independent domain 
approximation, one sees a clear discrepancy. 
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q 


Exact A(q) 


Exact B(q) 


Approximate A(q) 


Approximate B(q) 


1.2 


1.0769 


.153 


1.0675 


.133 


1.5 


1.1748 


.344 


1.1532 


.299 


2. 


1.3062 


.597 


1.2685 


.517 


3. 


1.4998 


.963 


1.4370 


.832 


5. 


1.7537 


1.439 


1.6630 


1.239 


10. 


2.1066 


2.101 


1.9777 


1.809 


100. 


3.2273 


4.369 


3.0065 


3.804 



for q < 2 and (|5jj52j) for q > 2 



Table 1: The exact values of A(q) and B(q) obtained from ( fls| , 
are compared to those coming from the independent domains approximation (28,29) 
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Second moment for various q's 
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< x > 




1.30 
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Logarithm of the integrated distribution of domains 




